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Supported magnetic nanoclusters: Softlanding of Pd clusters on a MgO surface 

M. Moseler 1 , H. Hakkinen 2 and U. Landman 2 
1 Theoretical Quantum Dynamics, Faculty of Physics, University of Freiburg, D-79104 Freiburg, Germany 
2 School of Physics, Georgia Institute of Technology, Atlanta, GA 30332-0^30 

(February 1, 2008) 

Low-energy deposition of neutral Pdjv clusters (N=2 — 7 and 13) on a MgO(OOl) surface F-center 
(FC) was studied by spin-density-functional molecular dynamics simulations. The incident clusters 
are steered by an attractive "funnel" created by the FC, resulting in adsorption of the cluster, 
with one of its atoms bonded atop of the FC. The deposited Pd2-Pd6 clusters retain their gas-phase 
structures, while for N>6 surface-commensurate isomers are energetically more favorable. Adsorbed 
clusters with N > 3 are found to remain magnetic at the surface. 

PACS: 68.47.Jn, 68.43.Bc, 75.70.-i 



Deposition of atomic clusters onto solid surfaces is a 
versatile surface-processing tool, with applications rang- 
ing from "micro-machining" and surface-smoothing to 
thin-film growth and fabrication of model nanocata- 
lysts Theoretical investigations, employing most 

often molecular dynamics (MD) simulations in conjunc- 
tion with semi-empirical classical interatomic potentials, 
guided many of the above experiments and provided 
valuable insights into the microscopic mechanisms of 
the deposition process p|-|lO|]. However, in cases where 
the dominant part of the cluster-surface interaction in- 
volves surface chemistry (that is, the creation or breaking 
of chemical bonds) [j7|, spin-dependent (magnetic) pro- 
cesses Q| , or surface defects of electronic origin (such as 
a F-center on an ionic surface) j3 11 Lq] , a full quantum 



description of the cluster deposition process is necessary. 

Here we report on a first-principles investigation of 
softlanding of Pdjv clusters (N=2-7 and 13) onto a 
MgO(OOl) surface containing a surface F-center, FC 
(oxygen vacancy) , yielding microscopic details about the 
cluster adsorption mechanism and the evolution of the 
atomic, electronic, and spin degrees of freedom during the 
process. We show that the interaction between the Pd 
cluster and the FC evolves from an initial long-ranged at- 
tractive polarization into chemical bonding involving the 
localized FC electronic state located in the MgO band 
gap. A common structural motif for the adsorption ge- 
ometry of the smaller clusters (A^<6) is found where the 
cluster retains its gas-phase geometry, whereas the larger 
clusters (N=7 and 13) adapt upon adsorption to the un- 
derlying surface rocksalt structure. The interaction with 
the surface quenches the spin for clusters with A^<3, re- 
tains the gas-phase triplet state (S=l) for 4 < N < 7, 
and for N=13 the gas-phase nonet (S=4) transforms to 
a septet (S=3) state. These results are of importance 
for understanding the activity of Pd/MgO nanocatalysts 
Jl3],[l4|], as well as for future investigations of supported 
magnetic clusters and nanoparticles. 

The Pd^/MgO system was treated within the frame- 



work of the local-spin-density functional (LSD) the- 
ory, with norm-conserving scalar-relativistic pseudo- 
potentials |ll| and self-consistent gradient corrections 
(PBE-GGA) |l| . For cluster impact energies below 
the band gap of the substrate material, the Born- 
Oppenheimer (BO) approximation provides a faithful de- 
scription of the collision dynamics and therefore we em- 
ployed the BO-LSD-MD method Q for the calculation 
of the electronic structure and the nuclear motion of the 
Pd cluster and the substrate. The MgO substrate with 
the F-center was modeled with a two-layer ab-initio clus- 
ter Mgi 3 0i2, embedded into a lattice of point-charges 
in order to simulate the long-range Madelung potential 
plU . The lattice parameter of the embedding part of 
the substrate was fixed to the experimental value (4.21 
A) of bulk MgO. The Pd^ cluster and the F-center's 
4 nearest-neighbor Mg atoms and 4 nearest-neighbor O 
atoms of the first layer were treated dynamically during 
the deposition. To model the heat conductivity of the 
extended MgO surface the equations of motion of the 
dynamic surface atoms included an added damping term 
with a damping constant ttujo/Q JTof , where u>d is the 
Debye frequency of bulk MgO. 

The initial spin states (triplet for N=2—7 and nonet 
for N=13) and geometrical structures of the Pd^ clus- 
ters were taken from our recent gas-phase study |lq ] . The 
clusters were placed with a random orientation 4 A above 
the FC (measured from the cluster atom closest to the 
surface) and an initial velocity directed perpendicular to 
the MgO surface, corresponding to a kinetic energy of 
0.1 eV per atom to simulate softlanding conditions. JToj ] 
The spin of the cluster-substrate system was dynamically 
evaluated at each MD time step. Subsequent to the dy- 
namical evaluation of the deposition process for about 1 
ps the simulation was stopped, and starting from the last 
recorded configuration a corresponding potential energy 
minimum was located by an energy-gradient optimiza- 
tion with variable spin; other spin-isomers (SPIs) were 
optimized (starting from the aforementioned optimal 
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FIG. 1. A Pels cluster impinges with 0.1 eV/atom kinetic energy on an F-center in a MgO (001) surface. Pd atoms are depicted 
as blue, Mg as green and O as red spheres, (a) Isosurfaces of the highest occupied up-spin molecular orbital (HOMO-1), the 
highest occupied down-spin orbital (HOMO) and the lowest unoccupied down-spin orbital (LUMO) of the initial configuration 
at t=0 (color coding distinguishes the sign of the wave function). Note, that the 75 up-spin and 73 down-spin orbitals are both 
occupied, (b) The corresponding local densities of states of the surface (blue area for up- and red area for the down-spin DOS) 
and of the cluster (red line for up- and blue line for down-spin DOS). The Fermi level, Ef = 0. We note that the DOS of the 
isolated surface and free cluster are essentially identical to that shown here, except for the first peak below Ef that corresponds 
to the long-range interaction discussed in the text, (c) Snapshots from the MD simulation recorded at the indicated times, (d) 
Time-evolution of the kinetic energy (red line) and the z-component of the CM velocity (blue line) of the cluster, (e) Evolution 
of the HOMO-LUMO gap (black line), of the eigenvalue energy difference £74^ — £751 (red line, see panel (a) for explanation 
of the orbital numbers), and of the z-component of the cluster's CM coordinate (blue dashed line). The triplet-singlet-triplet 
transition is indicated by the black arrows drawn on the time axis, (f) Isosurfaces of spin polarization density for the optimized 
triplet (S=l) and singlet (S=0) states. Yellow and purple denote excess of up and down spins, respectively, (g) Isosurfaces of 
the orbitals in (a) for the optimized adsorbed cluster, (h) The local DOS corresponding to the optimized cluster (color coding 
as in (b)). 



configuration), in order to explore the thermal stabil- 
ity of the lowest-energy SPI. 

The adsorption of a single Pd atom on top of the FC 



(£FC site) is characterized by a strong binding energy 
(3.31 eV) and a short equilibrium adsorption distance 
(1.65 A), compared to adsorption on-top of an oxygen 
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FIG. 2. Structural and magnetic size-evolution of sup- 
ported Pdjv clusters, (a) GS structures of Pdjv (N=2,3,4,6,7 
and 13). Color coding as in Fig. 1 except for Pdi3 where a sub- 
set of the Pd atoms is colored in yellow in order to highlight 
the Pd7 subunit (blue), (b) Size-evolution of the adhesion 
energy E a d (red filled squares), the binding energy per atom 
Eh for the supported (blue solid dots) and free (blue circles) 
clusters, and the HOMO-LUMO gap of the supported (green 
solid diamonds) and free (green open diamonds) clusters, (c) 
Size evolution of GS spin S (red diamonds) and the distance 
of the highest cluster atom to the surface (blue solid dots). 
The inset in (c) shows the SPI energies A_E (with reference 
to the GS configuration, S = 3 for the adsorbed cluster and 
S — 4 for the free one) and corresponding activation temper- 
atures T = 2AE/k/(3N — 6), of supported (green solid dots) 
and free Pdi3 (green circles) clusters. 



(tO) atom at the ideal MgO surface (1.16 eV and 2.17 



A). The bonding between the Pd atom and the FC in- 
volves the localized FC electronic orbital, located in the 
band gap of MgO (separated from the top of the valence 
band by 2.3 eV), and (mainly) the d(m = 0) orbital of 
the Pd atom. The attractive interaction to the F-center is 
rather long-ranged extending up to about 5 A above the 
surface; e.g., the interaction energy of a Pd atom placed 
5.2 A above the FC is 0.1 eV. This weak attraction is 
due to polarization of the d(m=0) valence orbital of Pd 
by the FC. Surprisingly, we found that none of the other 
adsorption sites for the Pd atom, lying in the vicinity of 
the F-center (e.g., on-top of the neighboring oxygen (fO), 
on-top of the neighboring Mg atom (tMg), the Mg-Mg 
bridge (MVlgMg), the Mg-0 bridge (6MgO) and the Mg- 
Mg-0 hollow site (ZiMgMgO)), are stable; i.e., optimiza- 
tion starting from any of these sites leads to a sponta- 
neous (barrierless) transition to the aforementioned tFC 
configuration. We conclude that the F-center acts as a 
rather wide attractive " funnel" for the Pd atom, extend- 
ing several A both laterally and vertically [^0| . This fun- 
neling effect steers the incident cluster and dominates the 
dynamics of the initial phases of the deposition process, 
as illustrated in the following for the representative case 
of a Pds cluster. 

When the Pds cluster is placed 4 A above the oxygen 
vacancy, the FC electronic state (located just below Ep) 
combines with d-orbitals of the closest Pd atom to form 
two bonding molecular orbitals (see the up-spin HOMO-1 
and the down-spin HOMO in Fig. la). All other orbitals 
(for example the lowest unoccupied orbital (LUMO) of 
Pds shown in Fig. la) remain to a large degree eigen- 
states of the separated systems, and consequently the 
corresponding density of states (DOS in Fig. lb) may be 
represented as a superposition of those of the bare surface 
and the gas-phase cluster. 

The long-range attraction between the cluster and the 
FC accelerates the lowermost Pd atom towards the tFC 
site (note the strong deformation of the Pds cluster at 
0.2 ps in Fig. lc and the increase of the kinetic energy 
in Fig. Id). Subsequently, other Pd atoms are attracted 
to neighboring &MgO positions (Fig. lc, t=0.5 ps) ac- 
companied by additional release of kinetic energy. Con- 
sequently, the center of mass (CM) velocity toward the 
surface increases to almost twice its initial value (Fig. Id) 
leading to a strong flattening of the cluster at t=0.5 ps 
(see the minimum in the z component of the cluster CM 
in Fig. le). The cluster shape deformation causes a tran- 
sient reordering of the molecular orbitals, i.e. it raises 
the energy of the up-spin HOMO-1 level (marked 75 1 in 
Fig. la) and turns it into a HOMO at t=0.32 ps, and 
even into a LUMO state at t=0.36 ps; this sequence is 
portrayed in Fig. le by closing of the HOMO-LUMO 
gap (black curve) and the minimum in the eigenvalue 
energy difference £751—6741 (red curve). Since the down- 
spin LUMO (74 J. in Fig. la) temporarily becomes the 
HOMO state, the total spin flips for a short time period 
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from S=l to S=0 (Fig. le). After 0.5 ps the cluster re- 
coils the reverse process drives the cluster back into the 
triplet spin state at t=0.63 ps (Fig. le). 

Optimization of the adsorbed cluster after a 1.2 ps 
MD simulation (see Fig. lc for the last MD configura- 
tion), resulted in a trigonal bipyramide structure (which 
coincides with the gas-phase optimal configuration) ly- 
ing with a triangular facet against a tFC- MvlgO- &MgO 
surface triangle (Fig. If). The spin polarization of the 
triplet ground state (see S — 1 isosurface in Fig. If) re- 
sembles that of the free cluster with a minor additional 
contribution from four surface oxygen atoms closest to 
the FC. As expected from our gas-phase calculations fl3] 
the slightly higher-lying singlet state (AE=2A meV) con- 
sists of an anti-ferromagnetic ordering of the local mag- 
netic moments (see 5 = isosurface in Fig. If). The 
spatial character of the orbitals close to Ep and the sur- 
face and cluster contributions to the DOS of the triplet 
ground-state of Pd5/MgO(FC) are shown in Figs, lg and 
lh, respectively. 

Using the above-mentioned methodology, we have de- 
termined the ground-states for the other deposited Pdjv 
clusters |2lj. For 3 < N < 6 we observed a regular 
structural size evolution (Fig. 2a) where the gas-phase 
GS structures are anchored to the MgO surface with 
one Pd atom on the £FC, another Pd on ZiMgMgO site 
(for N = 2) or 2 additional Pd atoms on MVlgO sites 
close to the tO position (for 3 < N < 6). However, 
for Pd7 and Pdi3 the free clusters transform to struc- 
tures that exhibit a higher degree of commensurability 
with the underlying surface, incorporating a Pdg and 
Pdy subunit, respectively (Fig. 2a). In this case, the 
loss in the intracluster cohesion is counterbalanced by 
a considerable gain of adhesion energy E a d [defined as 
E ad = J B(MgO(FC)) + J B(Pd 7V )- J B(Pd w /MgO(FC)), see 
the red curve in Fig. 2b]. Consequently, the cohesive en- 
ergy E c per Pd atom [defined as E c = (£(MgO(FC)) + 
NE(Pd) - £(PdAr/MgO(FC)))/A, see the blue curve in 
Fig. 2b] continues to increase after Pdg and remains well 
above the gas-phase E c values. 

The HOMO-LUMO gap (Fig. 2b) of the combined 
PdAr/MgO(FC) system is governed mainly by the metal 
cluster since the top part of its density of states lies in the 
MgO band gap. Most interestingly, the deposited Pd^r 
clusters with N > 4 remain magnetic: S=l for 4 < N < 7 
and S=3 for N=13. The crossover from nonmagnetic to 
magnetic states between Pd3 and Pd4 correlates with an 
increased "thickness" of the cluster (Fig. 2c), corroborat- 
ing our finding that flattening of the cluster on the surface 
tends to be accompanied by quenching of the spin (see 
discussion in the context of Fig. le). 

In general, the deposition of the cluster reduces the en- 
ergy separation between SPIs, thus lowering the thresh- 
old temperature for their coexistence. For instance, 
the triplet-singlet energy difference of supported Pd4 
is AE =65 meV compared to the gas-phase values of 



AE =136 meV; for Pdi3 five SPIs can be found within 
a 0.5 eV range, which expressed in terms of temperature 
corresponds to about 350 K (fig. 2c, inset). This result 
indicates that experiments aiming at distinguishing mag- 
netic states of the adsorbed clusters could be carried out 
at room temperature. 

In summary, we have used first-principles simulations 
in studies of the dynamics of soft-landing of Pd clusters 
on a MgO(lOO) surface containing an F-center defect, 
and for explorations of the structural and spin charac- 
teristics of the adsorbed clusters. The F-center creates 
an attractive "funnel" for the approaching metal clus- 
ter, resulting in preferred binding configurations with one 
Pd atom atop the F-center. A common structural motif 
was found for the smaller adsorbed clusters (Pd2 - Pd 6 ), 
where the cluster retains its gas-phase geometry. On the 
other hand, Pd7 and Pdi3 adapt upon adsorption to the 
underlying MgO rocksalt structure. Although the surface 
tends to reduce the spin of the adsorbed cluster, clusters 
larger than Pd3 remain magnetic at the surface, exhibit- 
ing several low-lying structural and spin isomers. These 
results provide the impetus for further investigations re- 
garding the interplay of structural and spin fluxionality 
of supported metal clusters and their catalytic properties, 
and call for an experimental verification of the predicted 
magnetic states of supported Pd clusters. 
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